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In this paper, we consider a spherical symmetric metric to extract the hydrostatic equilibrium 
equation of stars in (3 + 1)—dimensional gravity’s rainbow in the presence of cosmological constant. 

Then, we generalize the hydrostatic equilibrium equation to d-dimensions and obtain the hydrostatic 
equilibrium equation for this gravity. Also, we obtain the maximum mass of neutron star using the 
modern equations of state of neutron star matter derived from the microscopic calculations. It is 
notable that, in this paper, we consider the effects of rainbow functions on the diagrams related to 
the mass-central mass density (M-pc) relation and also the mass-radius (M-R) relation of neutron 
star. We also study the effects of rainbow functions on the other properties of neutron star such 
as the Schwarzschild radius, average density, strength of gravity and gravitational redshift. Then, 
we apply the cosmological constant to this theory to obtain the diagrams of M-pc (or M-R) and 
other properties of these stars. Next, we investigate the dynamical stability condition for these stars 
in gravity’s rainbow and show that these stars have dynamical stability. We also obtain a relation 
between mass of neutron stars and Planck mass. In addition, we compare obtained results of this 
theory with the observational data. 


I. INTRODUCTION 

In generalization from Galilean Relativity to Special Relativity, one has to modify the kinematic equations of 
Galilean Relativity to obtain an invariant velocity scale. The same method could be followed in order to derive a 
theory, containing an invariant energy scale which is Planck scale. Such modification is known as the doubly special 
relativity [l|. In other words, we take two upper limits for the particles probing a spacetime into account: the speed 
of light and the Planck energy. Due to modifications in the kinematic structure of the doubly special relativity, 
the energy-momentum conservation laws and the Lorentz symmetry group are altered. There are several theories 
suggesting that standard energy-momentum dispersion relation is modified in limit of the Planck scale which among 
them one can name; string theory Q , loop quantum gravity Q and non-commutative geometry . The generalization 
of the doubly special relativity to incorporate curvature is known as gravity’s rainbow [I|. 

The consideration of doubly special relativity, hence gravity’s rainbow, is related to the effects of quantum spacetime. 
In other words, it is proposed that the classical theory of the gravity is an emerging one from quantum degrees of 
freedom which in result causes the theory to be an effective one (a-Q- Such an effective theory has an energy 
dependent metric in order to describe different phenomena 0- The construction of the metric is determined by a 
particle probing it [0. 

On the other hand, the Einstein theory of gravity has fundamental problem in UV limit which indicates that it 
requires modification. It is not necessary to modify the action in order to have UV completion theory. Alteration 
of the metric describing the spacetime in some cases would be sufficient. Such an approach is employed in Horava- 
Lifshitz theory of the gravity as well. The gravity’s rainbow also enjoys this feature. In other words, it is an UV 
completion theory which in IR limit reduces to classical theory of the Einstein gravity. It is worthwhile to mention 
that by considering proper values for the energy functions of the gravity’s rainbow, this theory would yield a model 
of the Horava-Lifshitz which indicates that these theories are related. The UV completion theories require that their 
energy-momentum dispersion relation to be modified. Such modification is also observed in discrete spacetime [0, 
spacetime foam , spin-network in loop quantum gravity (LQG) Q , ghost condensation [0 and non-commutative 
geometry M . In addition, experimental observations also confirmed that such modification must be considered in the 
UV limit 0. 
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Several studies are conducted in context of black holes by consideration of the gravity’s rainbow |16l - l21] . It was 
shown that the thermodynamic properties of black holes are modified in gravity’s rainbow in such a way that it admits 
the existence of remnant for black holes. In evaporation of the black holes in gravity’s rainbow, the temperature of 
black holes goes to zero while its size is finite. Therefore, the black hole is not fully evaporated. Due to this feature, 
it is proposed that this model could provide a solution regarding the information paradox (2^ and the formation of 
naked singularity at the last stage of the black hole’s evaporation. The existence of remnant is observed for Kerr 
black holes, Kerr-Newman black holes in de Sitter space, charged AdS black holes, higher dimensional Kerr-AdS 
black holes and black saturn [2^. In addition, due to this phenomena, there exists an energy limit in which the 
mini black holes can be produced at LHC [2^. It is worthwhile to mention that the gravity’s rainbow also holds 
the usual uncertainty principle [2^ [2^. Recently, there has been an increasing interest in gravity’s rainbow 
The wormhole solutions in the context of gravity’s rainbow have been investigated in Ref. [32|. A study regarding 


the extended informal approach in 
quantum has been obtained in Ref. 


Jjravity’s rainbow is done in Ref. [33j| . In addition, the rainbow metric from 
(34| . On the other hand, the effects of g ravity’s rainbow on the thermodynamics 
of different models of black holes have been investigated in literature . Einstein theory of gravity is an effective 

theory which confronts specific problems in different regimes such as UV or in describing some phenomena such as 
increasing accelerating expansion of the Universe which was proven by the observation of high red-shift supernova 
[ 3 B - I 37 II and the measurement of angular fluctuations of cosmic microwave background fluctuations [ssl - l^ . One of 
the best candidates for modifications of Einstein theory is adding the cosmological constant to Einstein’s Lagrangian 

miii. 


The physics governing the structure of stars includes the interaction between gravitational force and internal pressure 
which results into an equilibrium state. This is known as hydrostatic equilibrium equation (HEE) which plays the 
essential role in studying the structure of stars. On the other hand, it was proven that in studying the compact objects 
such as the neutron and strange quark stars, consideration of curvature of space-time and generalization of Newtonian 
gravity to general relativity is necessary for describing different phenomena and making accurate predictions. The 
first attempt for obtaining Einsteinian HEE for stars was done by Tolman, Oppenheimer and Volkoff (TOV) J4j-l46l|. 
The TOV approach toward studying the structure of the stars has been employed by several authors pTI - l^ . In 
addition, the generalization to modified theories of gravity such as F{R), F{G) (Sdl - f^ and dilaton gravities are 
done (for more details see Refs, [bll - f^ b 

There has been an ongoing debate regarding the consideration of the higher dimensionality in studying different 
phenomena. There are several reasons which are supporting such generalization: First of all, in context of particles 
interactions, the necessity of such generalization is expressed for having consistent theories. On the other hand, 
it is essential to consider higher dimensionality in order to have an effective theory of superstring (67l - [^ . The 
phenomenology of higher dimensionality has proven to be richer [toI - I^ . As for the stars, a study regarding higher 
dimensional mass-radius relation for a star with the uniform density was done in Ref. [73|. The consideration of 
Kaluza-Klein model for describing the neutron stars was done in Ref. [tJ. This study highlighted the importance 
of higher dimensionality. In addition, the relativistic anisotropic stars were investigated, and the effects of higher 
dimensions we re p ointed out [^ . The TOV equation of Einstein-A gravity in higher dimensions was studied in 
the Refs. [t^ . l77l| . In this paper, we are interested in obtaining HEE for Einstein-A in gravity’s rainbow and its 
generalization to d-dimensional case. 

The outline of the paper will be as follows. First, we study a spherical symmetric metric and extract the HEE in 
Einstein-A gravity’s rainbow for (3 -I- 1)—dimensions. Next, we obtain a global equation of hydrostatic equilibrium 
for compact stars in the higher dimensions in this gravity. The maximum mass for neutron star in the presence and 
absence of cosmological constant in gravity’s rainbow is investigated. In the next section, we compare results obtained 
for theory with observational compact object. The last section is devoted to closing remarks. 


II. (3 -h 1) DIMENSIONAL HEE IN GRAVITY’S RAINBOW 


Here, we present the Einstein gravity with the cosmological constant. The action of this gravity is given by 

1 


Ig = - 


IGtt 


d‘^x^/^{R - 2A} -b / 


Matt 5 


( 1 ) 


M 


where R is the Ricci scalar, A is the cosmological constant and iMatt is the action of matter field. Varying the action 
0 with respect to the metric tensor the equation of motion for this gravity can be written as 


Gl + Kg'' = KTl 


( 2 ) 


where K = 


SttG 


and also, Gjj and are the Einstein and energy-momentum tensors, respectively. 


3 


Here, we intend to obtain the static solutions of Eq. Q in gravity’s rainbow. To do so, we use the metric which 
is energy dependent. As it was pointed out, one of the basics for gravity’s rainbow is deformation of the standard 
energy-momentum relation 


E^L^{e)-p^H^{e) =m^, (3) 

in which e = where Ep is the Planck energy. Considering the mentioned upper limit for energies that a particle 
can obtain, we have 


e<l. (4) 

The functions L{e) and H(e) are rainbow functions which satisfy the following conditions, 

lim L{e) = 1, lim i?(e) = 1. (5) 

e->0 E-s-O 

These conditions ensure the standard energy-momentum relation in the infrared limit. Now, we are in a position 
to construct the energy dependent metric as follows m 

h{£) (6) 


where 


eo(e) 



ei{e) 



(7) 


in which the tilde quantities refer to the energy independent frame fields. Therefore, one can assume a spherical 
symmetric space-time in the following form, 


ds^ 


fjr) 

L^e) 


dt^ 


1 

£f2(e) 


/ dr'^ 

Kdir) 


+ {d9^ + 



( 8 ) 


where /(r) and g{r) are radial dependent functions which should be determined. In addition, L(e) and i^(e) are the 
energy functions which will be introduced later. 

The energy-momentum tensor for a perfect fluid is 

= (cV -H P) - Pg^'', (9) 

where p and P are density and pressure of the fluid which are measured by local observer, respectively, and is 
defined as 




L{£) 




, 0 , 0,0 , 


which is the fluid four-velocity with following restriction 

= 1 . 


( 10 ) 


( 11 ) 


Using Eq. (0 and the metric introduced in Eq. ®, we can obtain the components of energy-momentum tensor 
for (3 -I- l)-dimensions as follows 


^0 — 


& Tl = = Ti = -P. 


( 12 ) 


We consider the metric ([5|) 
forms 


and Eq. m for perfect fluid and obtain the components of Eq. @ with the following 


Kc^r^p = Ar'^ + {1 — g) E[^{£) — rg'E[^{£), (13) 

Kr^fP = -Ar^f - {1 - g) H^{£)f + rgH^i£)f, (14) 

AKrfP = -AArf + 2 [gf) 'H^{£)f + r [g'f + 2gf''\ H^{£)f - rgH^{£)f'^ (15) 


where /, g, p and P are functions of r. We note that the prime and double prime denote the first and second 
derivatives with respect to r, respectively. 
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Using Eqs. (11311151) and after some calculations, we have 

dP {c^p + P) f 


dr 


2 / 


= 0 . 


Now, we obtain f from Eq. m as follow 


/' = 


[r^ (A + j^P) + (l-g) H^e)] f 
rgH'^{e) 


Then, for obtaining g, we use Eq. (US 


5=1 + 


A 


2 ATM, 


eff 




Attt 


(16) 


(17) 


(18) 


where Mgjf is a function of r and e: Mf.ff (r, e) = f It is notable that the obtained relation for effective 

mass is a function of radial coordinate and energy functions. This is a direct contribution of the gravity’s rainbow. By 
inserting Eqs. (HZl) and (fTSl) in Eq. we can extract the HEE in Einstein-A gravity’s rainbow for (3 + l)-dimension 
as 


dP _ [3c^GM^ffII^(e) + r3 (Ac^ + 127rGP)] (c^p + P) 
dr (?r [QGMeffH‘^{e) — (?r (Ar^ + 3i7^(e))] 


(19) 


We see that this equation is different from HEE, and for A = 0 and H(e) = 1, it yields the usual TOV equation 
Also, for A 7 ^ 0 and i7(£) = 1, it reduces to HEE obtained in Ref. fn\ . 

Considering different phenomenologies, one can employ three different cases for rainbow functions. 

Motivated by studies conducted in loop quantum gravity and non-commutative geometry, the first case has the 
following energy functions, Bizi 


L(e) = 1 & H{e) = — rye". 

Using Eqs. (|T^ and (l20l) . the HEE is 

dP [3c^GMe// [1 — rye”] + r^ (Ac'* + 127rGP)] (c^p + P) 


dr 


c^r [6GMe// [1 — rye"] — c^r (Ar^ + 3 [1 — rye”])] 


( 20 ) 


( 21 ) 


where for ?y = 0, Eq. (I2ip reduces to the HEE obtained for Einstein-A gravity in Ref. [^. Here, since H{e) is a 
square root function, one can find follo-wing condition regarding its parameters 


ry < 


1 


The second case is related to the hard spectra from gamma-ray bursts which has energy functions as 

- 1 


L{e) = 


/3e 


& H{e) = 1. 


( 22 ) 


It is notable that due to the structure of obtained TOV equation (Eq. (fT^ b this case of energy functions will 
lead to absence of the effects of gravity’s rainbow. In other words, for this case, the TOV equation is independent of 
energy functions which is not of our interest. 

The third case of energy functions is due to consideration of constancy of the velocity of light which leads to the 
following relation, [s^ 


L{e) = Hie) = 


1 — Ae 


(23) 


We use Eq. 


for obtaining the HEE, so we have 


dr 


3c^GMeff (Ac'* -h 127rGP) 


(c^p -I- P) 


-c2rfAr2 + 3(T^) 


( 24 ) 
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where for A = 0, Eq. (l24l) reduces to Einstein-A gravity [t^- Since the energy functions must be positive valued (for 
avoiding the change of metric signature), one can extract following condition for this model of gravity’s rainbow 


A 


1 


< 

£ 


III. HEE OF GRAVITY’S RAINBOW IN HIGHER DIMENSIONS 


Now, we are going to obtain the HEE in Einstein-A gravity’s rainbow for higher dimensions (d > 5 where d 
represents the dimensionally of space-time). We consider the following metric 


ds^ = ^ 


LHe) 


dr'^ 


\g{r) 


^dnt 


k ) 


where 


d-22-1 

dnl = d02 ^ JJ- 

i=2 j=l 

Then, we must obtain Eq. m for arbitrarily dimensions, therefore we have, 


(25) 


(26) 


T° = = r| = r| =... = = -p. 


(27) 


We use the above equation to obtain a global relation for the HEE in higher dimensions for Einstein-A gravity’s 
rainbow. Using Eqs. © and (EH) for the metric (1251) . the components of the equation of motion for Einstein-A 
gravity’s rainbow can be written as 

= Ar^ + ^d-2){d-2,) (28) 

K^fP = -Ar^f - id-^)id-^) (^ _ ^ , (29) 

^KarfP = -4Ar/2-^fc^M^^(l-5)ij2(£)/2 + ^[57' + 25/"]i7"(e)/ 

r 

+2 {d - 3) {gf) - rgH^{e)f^, (30) 

in which Kd = and Gd is the gravitational constant in d—dimensons and Gd is defined as Gd = GVd-i, where G 
denotes the four dimensional gravitational constant and also, Vd -4 is the volume of extra space. We obtain Eq. (HU) 
in the same manner as it was described for (3 -I- l)-dimensions, then using Eq. (1^ . we get /' as follows 


/' = 


(A -f KdP) + (1 _ g) H^e) 


rg{d — 2)H^{e) 


We calculate g of Eq. (1^ as 


_ 2A 2 c^dCdMeff {r,e)T {dL^) 

^ {d-l)id-2)H^ief (d-2)^(^-i)/2r^-3 


(31) 


(32) 


It should be noted that we have used Me/f (r, e) = f r{(d-i)/ 2 ')H‘^{e) ^p{r)dr in the above equation and also E is 
the gamma function, which satisfies some conditions as r(l/2) = r(l) = 1 and r(a; -I-1) = xT{x). 

Using Eqs. (131|) and (l32|) in Eq. (fTHll . we can get the HEE in the Einstein-A gravity for d-dimensions 


dr 


\{d-l){d-3)r(i^^)PKdM,ffH^{s) , 

(, {,d-l)KdP' 

)i 

(c^p -1- P) 


-1 ^ 

2 

r 

A 1 G-l) 

fr{^)c'^K^Meff (d-2)r‘^-3^ 

H^{e) 


^ J.d-1 


2 y 


( 33 ) 


where for (3 -I- 1)—dimensional limit, Eq. (1551) reduces to Eq. (1151) . 
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FIG. 1: Equation of state of neutron star matter (pressure, P (10^® erg/cm®) versus density, p (10^® gr/cm®)). 


IV. STRUCTURE PROPERTIES OF NEUTRON STAR 
A. Equation of state of neutron star matter 

It is stated that one can derive the properties of neutron star matter by obtaining the equation of state of neutron star 
matter. The constituents of the interior part of a neutron star include the neutrons, pro tons, electrons and muons 
which are in charge neutrality and beta equilibrium conditions (beta-stable matter) |8l| . On the other hand, for 
studying the equation of state of stars, there are several approaches which among them one can name the microscopic 
constrained variational calculations based on the cluster expansion of the energy functional [12,113 ■ approach, 

the new Argonne AV18 and charged dependent Reid-93 are employed as the two-nucleon potentials (see Refs. [84.[85l|- 
for more details). Also, a good convergence, lack of need for any free parameter in formalism and more accuracy 
comparing to other semi-empirical parabolic approximation methods are of the advantages of these methods. The 
modern nucleon-nucleon potentials are isospin projection (Tz) dependent in which a microscopic computation of 
asymmetry energy is carried on for the asymmetric nuclear matter calculations (see Ref. for more details). 

In this paper, we employ AV18 potential [H, 113 for calculating the modern equation of state for neutron star 
matter. Then, we will study some physical properties of neutron star structure. Our result for obtained equation 
of state of neutron star matter is presented in Fig. [TJ We extract the mathematical forms for the equation of state 
presented in Fig. [T]as 


P = Y^Ap'^-\ (34) 

where Ai are 

Ai = -3.518 X 10“®^ 

A 2 = 3.946 X 
As = -1.67x10"^®, 

A 4 = 3.242 X 10-1°, 

As = -1.458 X 10°, 

Ae = 2.911 X 10i°, 

A 7 = -9.983 X 10°i. 

In order to have a better insight regarding the properties of these neutron stars, we study both energy and stability 
conditions in the following subsections. 
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TABLE I: Energy conditions for neutron star. 



1 Bo(10®^; 

) NEC WEC SEC DEC 

3895.99 

1910.39 

/ / / / 



P 


FIG. 2: Sound speed (u^/c^) versns density (p x 10^® (gr/cm®)). 


1. Energy conditions 

Here, we investigate the energy conditions such as the null energy condition (NEC), weak energy condition (WEC), 
strong energy condition (SEC) and also, dominant energy condition (DEC) at the center of neutron star. We have 


NEC 

-)> 

Pc + Pc^ 0 , 


(35) 

WEC 

-)> 

Pc p Pc^ 0 , 

IV 

o 

(36) 

SEC 


Pc p PcP 0 , 

& iPc P PcP^, 

(37) 

DEC 


Pc P \Pc\ 1 


(38) 


where Pc and Pc are the density and pressure at the center of neutron star (r = 0), respectively. Considering Eig. [T] 
and the above conditions (|351I38L our results are given in table ID According to Eig. [T]and table HI we observe that all 
energy conditions are satisfied. So, the equation of state of neutron star matter introduced in this paper is suitable. 


2. Stability 

In order to investigate the stability of equation of state of neutron star matter for a physically acceptable model, 
one expects that the velocity of sound {v) be less than the light’s velocity (c) [H, By considering the stability 
condition (0 < < c^) and Eq. (|M1) . and comparing them with diagrams related to speed of sound-density 

relationship plot in Fig. [5J it is evident that this equation of state of neutron star matter satisfies the inequality 

0 < . 

Our investigations show that our equation of state of neutron star matter satisfies both energy and stability condi¬ 
tions. Now, we focus on investigation of gravitational mass and radius for neutron stars in gravity’s rainbow. 









TABLE II: Structure properties of neutron star for different values of H{e) (1.2 <H{e)< 1.67). 


Hie) Mmax {Mq) R {km) 

Rsch {km 

,) -p (10^®3 cm"®) 

a(lO-i) 

2(10"’^) 

1.70 

2.87 

14.32 

8.46 

0.46 

5.91 

5.63 

1.67 

2.81 

14.00 

8.28 

0.49 

5.91 

5.65 

1.60 

2.69 

13.47 

7.93 

0.52 

5.89 

5.59 

1.50 

2.52 

12.63 

7.43 

0.59 

5.88 

5.58 

1.40 

2.37 

11.79 

6.99 

0.69 

5.93 

5.66 

1.30 

2.19 

10.95 

6.46 

0.79 

5.90 

5.61 

1.20 

2.02 

10.10 

5.95 

0.93 

5.89 

5.61 


B. Properties of neutron stars in gravity’s rainbow without the cosmological constant 


Study of compact objects such as neutron stars has been of great interest for the astrophysicists in the last two 
decades (see Refs. [8^4^ for more details). Of the most important features of neutron stars are the maximum mass 
and radius, and determining them. Because, there is a critical maximum mass in which for masses smaller than that 
the degeneracy pressure originated from the nucleons prevents an object from becoming a black hole [Slj. Therefore, 
obtaining the maximum gravitational mass of neutron stars is of great importance in astrophysics. Due to many 
errors in direct ways of measuring mass and radius of the neutron stars by observations of the X-ray pulsars and 
X-ray bursts, one is not able to obtain an accurate mass and radius for these stars. On the other hand, using the 
binary radio pulsars [qgI - I^ . leads to highly accurate results for the mass and radius of neutron stars. 

Now, by employing the equation of state of neutron star matter which is presented in Fig. [1] and numerical 
approach for integrating the HEE obtained in the equation (11911 . we obtain the maximum mass and radius of neutron 
star. It is worthwhile to mention that the neutron star mass and radius are central mass density (pc) dependent in 
this approach. For this purpose, one can consider the boundary conditions P{r = 0) = Pc and m{r = 0) = 0, and 
integrates Eq. dUD outwards to a radius r = i? in which P vanishes for selecting a pc- This leads to the neutron star 
radius R and mass M = m{R). The results are presented in v arious tables and figures. 

The Einstein gravity case has been investigated in Ref. |lOCll| , and the maximum mass of neutron stars by using the 
modern equations of state of neutron star matter derived from microscopic calculations was obtained. It was shown 
that the maximum mass for neutron stars is about I. 68 M 0 in Einstein gravity. The effect of the cosmological constant 
on maximum mass of neutron star was investigated, and it was pointed out that considering the positive values of the 
cosmological constant, the maximum mass of neutron star decreases [za and also, the behavior of neutron star (for 
example; diagram related to the maximum mass-radius) for the negative values of the cosmological constant are not 
logical (see Ref. for more details). 

On the other hand, the maximum mass of neutron star is still an open question. There are some observational 
evide nces which indicate that the maximum mass of neutron stars can be more than 1.68Mq. For example, Jacoby 
et al |l0l| and Verbiest et al [I02j| for a binary system used the detection of Shapiro delay to me asur e the masses of 
both the neutron star and its binary companion. Also, using the same approach, Demorest et al |l03j | has performed 
radio timing observations for the binary millisecond pulsar PSR J1614-2230 where the measured mass for this pulsar 
was obtained abo ut 1.97Mq. The mass of othe r com pact objects were obtained about I^Mq for Vela X-I [lOJ, PSR 
J0348-f0432 [lo| about 2 .OIM 0 , 4U 1700-377 [lO^ about 2 . 4 M 0 , and J1748-2021B [lo3 about 2.7Mq. Also, A. W. 
Steiner et al determined an empirical dense matter equation of state from a heterogeneous dat a set of six neutron 
stars. They also showed that, the radius of a neutron star must be in the range R < (11 ~ 14) km [l08j| . In the present 
paper, we would like to see whether the gravity’s rainbow and obtained maximum mass for neutron stars by employing 
the modern equations of state of neutron star matter derived from microscopic calculations can predict existence of 
the maximum mass and radius of neutron stars more than I. 68 M 0 and less than 14fcm, respectively. Therefore, we 
consider the neutron stars with radius less than R < lAkm and separate the maximum mass for Mmax > 2 M 0 and 
J4t„ax < 2Mq in tables [TT] and IIIIl respectively. 

It is notable that, for Plje) = 1, the maximum mass of neutron star reduces to the result that was obtained in 
Einstein gravity [t^ . I100I| . as expected (because, this case (i7(e) = 1) is denoted as the absence of gravity’s rainbow. 


In other words, in this case, the effects of the gravity’s rainbow are vanished). On the other hand, our results show 
that, by increasing rainbow function more than 1 {Hie) > 1), the maximum mass of neutron star increases (see tables 
El and uni). In other words, our results cover the mass measurement of massive neutron stars, and also, predict that 
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TABLE III: Structure properties of neutron star for different values of H{e) (0.5 <H(e)< 1.1). 


H{£) Mmax (Mq) R (km) 

Rsch {km 

,) p (10^®3 cm-^) 

a(lO-i) 

a(lO-i) 

1.10 

1.85 

9.26 

5.45 

1.11 

5.88 

5.59 

1.00 

1.68 

8.42 

4.95 

1.34 

5.88 

5.58 

0.90 

1.51 

7.58 

4.45 

1.65 

5.87 

5.56 

0.80 

1.35 

6.73 

3.98 

2.10 

5.91 

5.64 

0.70 

1.18 

5.89 

3.48 

2.74 

5.91 

5.62 

0.60 

1.01 

5.05 

2.98 

3.72 

5.90 

5.61 

0.50 

0.84 

4.21 

2.48 

5.35 

5.89 

5.58 



R(km) 


FIG. 3: Gravitational mass versus central mass density (left diagram) and radius (right diagram), for H{e) = 0.9 (continuous 
line), H{e) = 1 (doted line), H{e) = 1.2 (dashed line), H{e) = 1.5 (dashed-dotted line) and H{e) = 1.7 (bold line). 


the mass of neutron stars in gravity’s rainbow can be in the range larger than I. 68 M 0 (see tables HIl and IIIII for more 
details). In order to conduct further investigations, we plot diagram related to the gravitational mass versus central 
mass density (or radius) in Fig. [3l The variation of the gravitational mass versus different values of the rainbow 
function is presented in this figure. On the contrary, for H{e) < 1, both maximum mass of neutron stars and radius 
of these stars are smaller than the Einstein case (see table Hill and Fig. [3]). This emphasizes the contributions of the 
gravity’s rainbow on properties of neutron stars. 

In the above tables, we used the obtained TOV equations in Eqs. (1211) and (IMll for H{e) < 1 and H{e) > 1, 
respectively, to extract the structure of neutron stars. As one can see, due to the structure of obtained TOV equation 
(Eq. (fT31) ). for case (I33|) . the obtained TOV equation is independent of energy functions which is not of our interest. 

In the following, we are going to investigate other properties of neutron star in the gravity’s rainbow. 


1. Schwarzschild radius 

Using Eq. mi), and this fact that g{r = Rsch) = 0, we calculate Schwarzschild radius for the obtained masses in 
gravity’s rainbow as 


RSch 




(39) 


As it was pointed out, the mass of neutron star depends on rainbow function and thus we expect the Schwarzschild 
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radius varies depending on choices of rainbow function. The results are presented in tables |TT] and IIIII It is evident 
that, the Schwarzschild radius is an increasing function of rainbow function. 


8. Average density 

The average density of a neutron star has the following form 

iMeff 


(40) 


where it is a decreasing function of H{e). 

Obtained central density by using the mentioned equation of state of neutron star matter is about 3.9 x 10^® g cm~^ 
and this density is larger than the normal nuclear density, po = 2.7 x 10^^g cm~^ |l09l |. On the other hand, by 
considering the rainbow function less than 0.6 {H{e) < 0.6), obtained average density for neutron stars is larger 
than the central density (see the last row of table Hill) . In other words, when the rainbow function is larger than 0.6 
{H{e) > 0.6), the neutron star is compatible. 


3. Compactness 

The compactness of a spherical object may be defined as the ratio of Schwarzschild radius to the radius of object 
(a — Rsch / R), which may be indicated as the strength of gravity. Obtained a in this gravity shows that the strength 
of gravity is almost the same for these neutron stars. 


4- Gravitational redshift 


Using Eq. dn, we can obtain the gravitational redshift in the following form 


z = 


1 



2GMeff 

c^R 


- 1 , 


(41) 


where M^ff is the mass of neutron stars in which is rainbow function dependant. This relation indicates that to have 
a physical redshift, the radius of a neutron star should be greater than the Schwarzschild radius, R > = Rg^j^. 

The results related to the gravitational redshift shows that for neutron stars with different mass, this quantity is 
almost the same (see the last column of tables m and mi. We plot the diagrams related to the gravitational redshift 
versus mass in Fig. 0] This figure shows that, for different values of rainbow function, this quantity is not affected 
considerably. 


5. Dynamical stability 


The dynamical st abili ty of the stellar model against the infinitesimal radial adiabatic perturbation was introduced 
by Cha ndrasekha r [ng. Then, this stability condition was developed and applied to astrophysical cases by many 
authors [UMii. The adiabatic index (7) is defined in following form 


pc^ + P dP 
c^P dp 


(42) 


In order to investigate the dynamical stability condition, 7 should be more than f (7 > f = 1.33) everywhere 
within the isotropic star. For this purpose, we plot two diagrams related to 7 versus radius for different values of 
rainbow functions in Figs. [5]and|6l As one can see, this stellar model is stable against the radial adiabatic infinitesimal 
perturbations. 

Also, we plot the pressure (density) versus distance from the center of neutron star. As one can see, the pressure 
and density are maximum at the center and decreases monotonically towards the boundary (see Figs [7] and [8|. 

As o ne c an see, considering the radius of neutron stars less than lAkm (according to the work of A. W. Steiner 
et al. [l08j| l and by using the average density, we obtained an upper and a lower limitations for rainbow function. 
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FIG. 4: The gravitational redshift versus mass for H{e) = 1.1 (continuous line), H{e) = 1.3 (doted line) and H{e) = 1.5 
(dashed line). 



FIG. 5: Adiabatic index versus radius for H{e) = 1.80 (dashed line), H{e) = 1.60 (continuous line) and H{e) — 1.4 (dashed- 
dotted line). 


Therefore, the obtained solutions in gravity’s rainbow are valid in the range of 0.6 < H{e) < 1.67. In other words, 
gravity’s rainbow cover all neutron stars and pulsars which have the maximum mass less than 2 . 8 IM 0 (see tables HIl 
and ED). Indeed, these results cover wide range of neutron stars and pulsars such as Vela X-1, PSR J0348-I-0432, 4U 
1700-377, J1748-2021B, PSR J1903 -h 327, Gen X - 3, PSR B1913 -h 16, PSR J0737 - 3039, PSR J0737 - 3039B and 
SMC X - 1. 


C. Relation between the mass of neutron stars and Planck mass 


Here we derive the relation betw een neutron star mass and the Planck mass in the gravity’s rainbow according to 
Burrows and ©striker’s work [ll5j| . We know that in addition to the degenerate pressure of nucleons, the pressure 
due to the strong repulsive inter-nucleons force balances the pressure due to the gravitation force in the neutron stars. 
Here, the nucleon-nucleon interaction which is so strong, is taken place through the pion exchange. Using this idea, 
we can consider the average density of a neutron star in term of the nucleus density using the following relation 



Pnuc 


(43) 







12 



FIG. 6: Adiabatic index versus radius for H{e) = 1.1 (dashed line), H{e) = 1.0 (continuous line) and H{e) = 0.9 (dashed-dotted 
line). 



FIG. 7: Pressure (erg/cm®) versus radius for H{e) = 1.80 (dashed line), H{e) = 1.60 (continuous line) and H{e) = 1.4 
(dashed-dotted line). 



FIG. 8: Density {gr/cm^) versus radius for H{e) = 1.80 (dashed line), 77(e) = 1.60 (continuous line) and 77(e) = 1.4 (dashed- 
dotted line). 
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where rup is the proton mass and is the Compton wavelength. Here is the pion mass. On the other 

hand, the mass of a ne utron star is set by a general relativity instability. The condition for Einstein gravity instability 
is R = (see Ref. |ll5l | for more details). Now, we use this idea for gravity’s rainbow. The Schwarzschild radius 
and the average density of a neutron star in gravity’s rainbow are obtained by equations (15^ and (HOI) as follows 


where 


RSch — 


P = 


Meff (r, e) 


2GMeff 2GM 

(44) 

•iMgf f 3M 

47ri?3 4:TTH^(e)R^^ 

(45) 

= ^2(£) / 47rr2p(r)dr. 

(46) 


Therefore, we have M^f / = H^e) ■ Now, we can derive the corresponding mass using the equations and (H51) and 
considering pnuc as follows 


M 



(47) 


where Mch is the Chandrasekhar mass. In Eq. (ITTl) . Met i^piVpi Vp = and where 

mpi is the Planck mass (see Ref. [Hi for more details). Therefore, the relation between the neutron star mass and 
the P lanck mass is obtained by above equation. It is notable that for i7(e) = I, Eq. (1471) reduces to Einstein gravity 
[HI. On the other hand, we can rewrite the equation (I47|) in the following form 

M ^ H^{e)-^EN. (48) 


where ^en = mpirjp 



D. Properties of neutron stars in gravity’s rainbow with cosmological constant 

In this subsection, we take the effects of variation of the cosmological constant into account. In other words, the 
effects of variation of the cosmological constant on maximum mass and radius of the neutron stars in the presence of 
Einstein gravity’s rainbow are investigated. Our results are presented in table IIVI 

As one can see, interestingly, the maximum mass of neutron star is a decreasing function of the cosmological 
constant (see table IIVl for more details). For further investigations of the behavior of the mass versus radius, we have 
plotted the following diagrams (Figs. and ITUl). 

It is evident from Figs. [5] and (TU] that the behavior of mass as a function of the central mass density (or radius) 
is highly sensitive to the variation of cosmological constant. For large values of cosmological constat (Fig. ITOl) . the 
behavior is completely modified. In this case, in the presence of large values of cosmological constant, the behavior 
of the mass versus radius of neutron star will be modified into a quark like behavior. 

It is a matter of calculation to show that the Schwarzschild radius of gravity’s rainbow in the presence of the 
cosmological constant is 


HHe) 


( 3GM,ff 


HHe) 




RSch — 


1/3 


Ai/3 
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SGMefj 

-^3- 


77^/3 (g) 



9G-Mipf\ 
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TABLE IV: Structure properties of neutron star for H{e) = 1 (up table) and H{e) — 1.5 (down table), respectively, with 
different values of A. 


A 

(m-^: 

) 

Mmax (Mq) R {km) 

Rsch {km 

') P cm"®) 

a(lO-i) 

z(lO-i) 

1.00 

X 

10" 

-16 

1.68 

8.42 

4.95 

1.34 

5.88 

5.58 

1.00 

X 

10" 

-14 

1.68 

8.42 

4.95 

1.34 

5.88 

5.58 

5.00 

X 

10" 

-14 

1.68 

8.41 

4.95 

1.33 

5.88 

5.58 

1.00 

X 

10" 

-13 

1.67 

8.40 

4.92 

1.33 

5.86 

5.54 

5.00 

X 

10" 

-13 

1.62 

8.34 

4.77 

1.32 

5.72 

5.29 

1.00 

X 

10" 

-12 

1.56 

8.25 

4.60 

1.31 

5.57 

5.03 

5.00 

X 

10" 

-12 

1.12 

7.47 

3.30 

1.27 

4.42 

3.38 

1.00 

X 

10" 

-11 

0.78 

6.65 

2.29 

1.25 

3.46 

2.36 


A (m ^) Mmax (Mq) R (km) Rsch {km) p (10^'^ff cm ^) a-(10 ^) z{10 ^) 


1.00 X 10"^® 

2.52 

12.63 

7.43 

5.94 

5.88 

5.58 

1.00 X 10"^'‘ 

2.52 

12.63 

7.43 

5.94 

5.88 

5.58 

5.00 X 10"^'‘ 

2.51 

12.62 

7.40 

5.94 

5.86 

5.54 

1.00 X 10"^® 

2.50 

12.61 

7.37 

5.92 

5.85 

5.51 

5.00 X 10"^® 

2.44 

12.51 

7.19 

5.91 

5.75 

5.34 

1.00 X 10"^® 

2.34 

12.38 

6.90 

5.86 

5.57 

5.03 

5.00 X 10"^® 

1.68 

11.21 

4.95 

5.69 

4.41 

3.38 

1.00 X 10"“ 

1.17 

9.98 

3.45 

5.58 

3.46 

2.36 


TABLE V: Properties of neutron star with and without the cosmological constant. 


H{e) 

A (m ^) 

Mmax {Mq) R {km) 

Rsch {km 

) P (10^®5 cm"®) 

a(10"i) 

z(10"i) 

1.00 

1 X 10"“ 

1.68 

8.42 

4.95 

1.34 

5.88 

5.58 


0 

1.68 

8.42 

4.95 

1.34 

5.88 

5.58 

1.30 

1 X 10"®® 

2.19 

10.95 

6.46 

0.79 

5.90 

5.61 


0 

2.19 

10.95 

6.46 

0.79 

5.90 

5.61 

1.60 

1 X 10"®^ 

2.69 

13.47 

7.93 

0.52 

5.89 

5.59 


0 

2.69 

13.47 

7.93 

0.52 

5.89 

5.59 


The above equation shows that the Schwarzschild radius is modified in the gravity’s rainbow. In other words, 
the Schwarzschild radius depends on the rainbow function. It is notable that, when H{e) = 1, the Schwarzschild 
radius reduces to the obtained Schwarzschild radius in Ref. The results show that the Schwarzschild radius is a 
decreasing function of the cosmological constant (see table lIVl for more details). We present the results of the average 
density and gravity strength in the presence of the cosmological constant in table IIVI We find that increasing the 
cosmological constant leads to decreasing both the strength of gravity and the maximum mass of neutron star. On 
the other hand, calculations related to the average density show that by increasing A, the average density decreases. 

We extract the gravitational redshift in gravity’s rainbow and in the presence of cosmological constant in the 
following form 


z = 


1 





2GM^ff 

c‘-‘R 


- 1 . 


(49) 


The results show that, the gravitational redshift is a decreasing function of the cosmological constant and an 
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FIG. 9: Gravitational mass of neutron star versus central mass density for H{e) 
A = 1 X 10“^® (doted line) and A = 1 x 10“^^ (dashed line). 


1.5, A = 1 X lO"^"* 


(continuous line), 



FIG. 10: Gravitational mass of neutron star versus radius for H{e) = 1.5, A = 1 x 10 A = 1 x 10 and A = 1 x 10 
from left to right, respectively. 


increasing function of Me//. In addition, in order to have a finite and real redshift, we should set 1+ > 0 

and therefore, R > Rsch, as we expected. 

Also, the dynamical stability for neutron stars in gravity’s rainbow and in the presence of cosmological constant 
shows that these stars have dynamical stability in all over the neutron stars (see Figs. [TT]and[T21 for more details). 

In order to investigate internal structure of the neutron star in more details, we plot the pressure (density) versus 
distance from the center of neutron star in the presence the cosmological constant in gravity’s rainbow. Figures [T51 and 
[14] show that, the pressure and density are maximum at the center and decrease monotonically towards the boundary. 

The value of the cosmological constant is an open question. On the other hand, from the perspective of cosmology, 
this value is about 10“^^ As one can see, the cosmological constant has no significant effect when its value 

is about 10“^^ m~‘^ (see table El for more details). In order to examine the effects of the cosmological constant on 
properties of the neutron star, we should consider a toy model in which its value is about less than 10“^"^ m“^. The 
results show that by decreasing the cosmological constant (less than 10“^^ m“^, A < 10“^'^ m“^), the maximum 
mass and also the radius of this star are not modified. In other words, for A < 10“^^ m“^, the cosmological constant 
does not affect the maximum and radius of neutron star (see table IIVI for more details). One of the results of this 
work is that when the value of the cosmological constant is about 10“^^ this constant does not play a role in 
the structure of neutron stars, but by taking larger values for it, the maximum mass and its radius are reduced. 
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•■(cm) 


FIG. 11: Adiabatic index versus radius for A = 1 x 10 H{e) = 1.80 (dashed line), H(£) = 1.60 (continuous line) and 
H{e) = 1.40 (dashed-dotted line). 



FIG. 12: Adiabatic index versus radius for H{£) = 1.40, A = 1 x 10 (dashed line), A = 5 x 10 (continuous line) and 
A = 1 X lO”*^"^ (dashed-dotted line). 



•■(cm) xio'" 

FIG. 13: Pressure (erg/cm^) versus radius for A = 1 x 10“^^, H{e) = 1.80 (dashed line), H{e) = 1.60 (continuous line) and 
H{£) = 1.40 (dashed-dotted line). 
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FIG. 14: Density (gr/cm®) versus radius for A = 1 x 10 , H{e) = 1.80 (dashed line), H{e) = 1.60 (continuous line) and 

H{e) = 1.40 (dashed-dotted line). 


TABLE VI: Properties of neutron star with and without the cosmological constant. 


H{e) Dimensions 

Mmax (Af©) 

R {km 

) Rsch (km) 

p (10^® gr cm ®) 

a 

1.70 

5 

9.22 X 10® 
32.06 X 10^^ 

17.52 

107.41 

0.39 

6.13 


6 

20.58 

341.98 

0.33 

16.62 

1.40 

5 

6.26 X 10® 
17.91 X 10^'^ 

14.43 

88.50 

0.58 

6.13 


6 

16.95 

281.65 

0.48 

15.72 

1.10 

5 

3.86 X 10® 

8.68 X lO^’^ 

11.33 

69.09 

0.94 

6.10 


6 

13.32 

221.23 

0.78 

16.73 

0.80 

5 

2.04 X 10® 
3.34 X 10^^ 

8.24 

50.52 

1.78 

6.13 


6 

9.68 

160.91 

1.48 

16.62 


V. HIGHER DIMENSIONS 

In this section, we extend our solutions to 5 and 6 dimensions to obtain the maximum mass of neutron stars by 
using the mentioned equation of state of neutron star matter. 

First, we consider the TOV equation for gravity’s rainbow without the cosmological constant. The results of our 
calculations show that the maximum mass of the neutron stars in 5 and 6 dimensions is not reasonable. In other 
words, the obtained maximum mass of neutron star is about 10 ®Mq and IO^^Mq for 5 and 6 dimensions, respectively 
(see table IVII for more details). In order to get a better picture of the results, we calculate the Schwarzschild radius 
for the obtained masses in gravity’s rainbow at these dimensions. The Schwarzschild radius for this gravity in higher 
dimensions is 


Rsch 


C^{d-2)TT^ ) 


(50) 


Also, we discuss the average density of the neutron star in 5 and 6 dimensions, in which can be written as 


P = 


2 Me 


15M 


t 


eff 




d = 5 
d = 6 


(51) 


The results are collected in table ED These results show that for 5 and 6 dimensions, obtained radius of the neutron 
stars in gravity’s rainbow are within the Schwarzschild radius. In other words, these objects can not be neutron stars. 
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TABLE VII: Mass and radius of neutron star through the observations and theory. 


name 

Observations 
M (Mo) 

R (km) 

H{e) M 

Theory 

max (-^^©) 

R {km) 

PSR J0348 -t 0432 

2.01 

13(±2) 

1.195 

2.01 

10.06(±0.01) 

PSR J1614 - 2230 

1.97 

13(±2) 

1.172 

1.97 

9.88(±0.01) 

417 1608 - 52 

1.74 

9.3(±1) 

1.035 

1.74 

8.71(±0.01) 


TABLE VIII: Predicted radius for neutron star through the theory. 



Observations 



Theory 


name 

M (Mo) 

R (km) 

H{e) 

M (Mo) 

R (km) 

J1748 - 2021R 

2.70 

unknown 

1.605 

2.70 

13.52(±0.01) 

4U 1700 - 377 

2.40 

unknown 

1.425 

2.40 

12.00(±0.01) 

PSR J1903 -t 327 

1.67 

unknown 

0.994 

1.67 

8.37(±0.01) 

Cen X — 3 

1.49 

unknown 

0.888 

1.49 

7.48(±0.01) 

PSR R1913-t 16 

1.44 

unknown 

0.856 

1.44 

7.20(±0.01) 

PSR J0737 - 3039 

1.35 

unknown 

0.802 

1.35 

6.75(±0.01) 

PSR J0737 - 3039R 

1.24 

unknown 

0.738 

1.24 

6.21(±0.01) 

SMC X - 1 

1.04 

unknown 

0.738 

1.24 

6.21(±0.01) 


but black holes. So the mentioned equation of state of neutron star matter could not anticipates the existence of 
neutron stars in higher dimensions. Therefore, the modern equations of state of neutron star matter derived from 
the microscopic calculations is a valid equation of state in 4 dimensions for neutron stars, but for higher dimensions, 
this equation of state should be modified. On the other hand, B. C. Paul et al showed that for anisotropic stars, 
by considering higher dimensions, the maximum mass of compact objects increases when dimensions of spacetime 
increase in which however attains a maximum in 5—dimensions, thereafter it decreases as one increases number of 
extra dimensions (see [t^, for more details). As it was pointed out before, these different results may be due to 
different equation of states or various kinds of compact objects (anisotropic and isotropic). So it will be interesting 
to study the equation of state of neutron star matter in higher dimensions. 


VI. THEORY AND OBSERVATIONS 

Here, we compare our results of gravity’s rainbow with observational data. For this purpose, we present these 
results in table IVIII 

As one can see, the results extracted in this theory match with the results obtained through the observations, with 
the exception that the radius obtained through gravity’s rainbow are smaller than the radius of compact objects. This 
is due to the fact that we have considered static neutron stars in our studies, so the radius of neutron stars is smaller 
than the radius of observational compact objects. 

Evidently for static cases of observational objects, our theory under consideration predicts a set of radii which 
for cases such as J1748-2021B and PSR J0348+0432, predicted radii are more than 13.52 (km) and 10.06 (km), 
respectively (see table IVIlTl and Fig. [T5]for more details). 


VII. CONCLUSIONS 

In this paper, by taking into account three models of gravity’s rainbow with (3+l)-dimensional spherically symmetric 
space-time, we studied the hydrostatic equilibrium equation (HEE). Then, we generalized our solutions to arbitrary d- 
dimensions (d ^ 5). Next, by employing obtained HEE, we conducted a study regarding the structure of neutron stars. 
Regarding microscopic calculations, we used the modern equations of state of neutron star to obtain the properties 
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R{km) 


FIG. 15: Gravitational mass versus central mass density (left diagram) and radius (right diagram), for PSR J0737 — 3039B 
(continuous line), PSJi J0737 — 3039 (dashed line), PSP i?1913 + 16 (doted line), 41/ 1608 — 52 (bold-continuous line), 
PSP J1614 — 2230 (bold-dashed line) and PSP J0348 — 0432 (bold-doted line). 


of neutron star. We have studied the effects of rainbow functions on the diagrams related to the mass-central mass 
density relation and also the mass-radius relation of neutron star. We showed that, for i?(e) > 1, by increasing H{e), 
the maximum mass of neutron star increases to more than 2Mq and for H{e) < 1 , the maximum mass is a decreasing 
function of the H{£). 

Then, we added the cosmological constant into this gravity and examined the diagrams related to M-p^ and M-R. 
It was pointed out that for A < 10“^'^ the cosmological constant (A) does not affect the maximum mass and 
radius of neutron stars. On the other hand, for A > 10“^^ m“^, the mass versus radius of a neutron star may present 
a quark-like star behavior (see Fig. [101 for more details). 

Next, we examined the effects of rainbow functions and the cosmological constant on the other properties of neutron 
star such as; the Schwarzschild radius, average density and strength of gravity. Also, we found a bound for rainbow 
function, so that for the value less than 0.6 {H{e) < 0 . 6 ), obtained average density for neutron stars was larger than of 
the central density, therefore, valid values for rainbow function in case of neutron stars are those that are larger than 
0.6 {H{e) > 0.6). Another interesting result of this paper is the effects of the cosmological constant on the so-called 
gravity strength. We found that for positive of the cosmological constant, the strength of gravity is a decreasing 
function of the cosmological constant which leads to decreasing the maximum mass of neutron stars. 

Also, we have obtained the gravitational redshift and showed that it was modified in this gravity. The gravitational 
redshift of neutron stars in this gravity was in the range z < 0.56 6 which is lower than the upper bound on the surface 
redshift for a subluminal equation of state, i.e z = 0.851 [ll6j |. In addition, we have investigated the dynamical 
stability and found that these stars were stable against the radial adiabatic infinitesimal perturbations. Our results 
showed that, the neutron stars with a mass more than 2Mq, can exist in gravity’s rainbow. Then, we compared our 
results with the empirical data and saw that our results were compatible with those of observations. Finally, using 
gravity’s rainbow, we obtained the radius of observational compact objects. We concluded that it is possible to obtain 
the mass and radius of compact objects by this theory through fine tuning. 

One of the most important results of this paper is regarding determining and finding an upper and a lower limitations 
for rainbow function. In other words, gravity’s rainbow was valid in the range of 0.6 < H(e) < 1.67. Furthermore, 
this gravity predicts a maximum limitation for the mass of neutron star which was about 2.81Mq. Indeed, gravity’s 
rainbow covers all neutron stars and pulsars which have the maximum mass less than 2 . 8 IM 0 . 

Using gravity’s rainbow we obtained relation between the maximum mass of neutron stars and the Planck mass. 
We showed that, this relation is modified and it depends on rainbow function. It is notable that, in the absence of 
rainbow function {H{e) = 1), this relation reduces to Einstein gravity. 

Considering the gravity’s rainbow and the obtained modified TOV in this gravity, we can investigate other compact 
objects such as white dwarf and quark stars. In addition, it is worth studying the effects of other equation of states 
on structure of the compact objects. On the other hand, anisotropic compact objects |ll7l - ll24 | in the context of 
gravity’s rainbow are also interesting subjects. We leave these issues for future works. 
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